The first quadrupole excitations in spherical nuclei and nuclear pairing 
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Abstract. Excitation energies and transition probabilities of the first 2 + excitations in even lead, tin and nickel 
isotopes are calculated within the self-consistent Theory of Finite Fermi Systems based on the Energy Density 
Functional by Fayans et al. A reasonable agreement with available experimental data is obtained. The effect of 
the density dependence of the effective pairing interaction is analyzed in detail by comparing results obtained 
with volume and surface pairing. The effect is found to be noticeable, especially for the 2 + -energies which are 
systematically higher at 200-300 keV for the volume paring as compared with the surface pairing case, the latter 
being in a better agreement with the data. 
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1 Introduction 

Since the famous paper of 1960 by Belyaev |[T1, it became 
clear that the low-energy nuclear dynamics is governed 
mainly with an interplay of pairing and the quadrupole 
degree of freedom, properties of the quadrupole excita- 
tions themselves depending crucially on the pairing. It is 
clearly seen directly from the experimental data. Indeed, 
there is no low-lying collective 2 + -state in magic nuclei 
where pairing is absent. E.g., in the 208 Pb nucleus the first 
2 + -state has the excitation energy ^2=4.08 MeV, but in the 
neighboring isotope 206 Pb it decreases to ^2=0.80 MeV. It 
is better to say that the low-lying 2 + -states and pairing have 
the common grounds, the unclosed shell structure where 
low energy virtual transitions of positive parity dominate, 
and pairing is just an approximate way to take into account 
the most important part of these correlations. In the paper 
0, we put attention to effects of the density dependence 
of pairing on characteristics of 2 + -states. 

Almost 50 years ago, in Ref. 0, the density-dependent 
ansatz for the effective pairing interaction was proposed 
within the theory of finite Fermi systems (TFFS) |4| with 
two (dimensionless) parameters y; n and y ex which deter- 
mine the strength of pairing interaction inside and out- 
side nucleus correspondingly. Arguments in favor of the 
version with the surface dominance, |yi n /y e xl - 10, were 
found in and later in [5] on the base of the detailed 
analysis of double odd-even mass differences in spherical 
nuclei. For brevity, we name it the "surface pairing", in 
contrast with the "volume pairing" model corresponding 
to equality y m -y ex . Very strong support of the surface ver- 
sion was produced by Fayans et al. |6| with the analysis, 
within the Energy Density Functional (EDF) method, of 
the odd-even staggering effect in charge radii of many iso- 
topic chains. The surface pairing version is supported also 
with the microscopic theory of pairing [7| starting from a 
realistic AW potential. Here, following to (2), we demon- 



strate that the systematic, partially new, data on the exci- 
tation energies a>2 in semi-magic nuclei evidence also in 
favor of the surface pairing. 



2 Brief formalism 

We describe the quadrupole excitations in spherical nuclei 
within the self-consistent TFFS [ 8 1 with the use of the EDF 
method of 1 6 1 . In this method, the ground state energy of a 
nucleus is considered as a functional of normal and anoma- 
lous densities, 



J S[p„(r), Pp (r),v„(r),v p (r)]d 3 



(1) 



This approach is a generalization for superfluid finite sys- 
tems of the original Kohn-Sham EDF method 1 9 1 . The nor- 
mal part of the EDF fi norm contains the finite-range central, 
spin-orbit and effective tensor nuclear terms and Coulomb 
interaction term for protons. We use the DF3-a functional 
IflOl which is a slight modification of the DF3 functional 
of ||6| concerning the spin-orbit and effective tensor terms. 

Within the TFFS, the response of a nucleus to the ex- 
ternal quadrupole field Vo exp (io)t) can be found in terms 
of the effective field. In systems with pairing correlations, 
equation for the effective field can be written in a compact 
form as 

V(w) = V (w) + rA(w)V(w), (2) 

where all the terms are matrices. In the standard TFFS no- 
tation 0, we have: 



V = 



r = 
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JXcS) MiM M 2 ((o) ) 
0(a)) -NiipS) N 2 (oj) 
O(-af) ~Ni(-<o) N 2 (-a>), 



(5) 



where £, Mi, and so on stand for integrals over e of the 
products of different combinations of two Green functions 
G(s) and two Gor'kov functios F^ l \e) and F (2) (e). They 
can be found in |4| and we write down here only the first 
of them which is of the main importance for us, 

/d 
— [G(s)G(e + a))- F (r> (s)F (2 \s + a))] . (6) 

Isotopic indices in Eqs. (f3]|6]l are omitted for brevity. In 
Eq. (O, T is the usual Landau-Migdal amplitude, 



Sp 2 ' 

is the effective pairing interaction, 

6 2 & 



5v 2 ' 



(7) 



(8) 



whereas the amplitudes = stand for the mixed 
second derivatives, 



6 2 B 
dpSv 



(9) 



In the case of volume pairing, we have - 0. The 
explicit form of the above equations is written down for 
the case of the electric (f-even) symmetry we deal with. A 
static moment of an odd nucleus can be found in terms of 
the diagonal matrix element (Ao\V(a) = 0)\Aq) of the effec- 
tive field over the state Aq of the odd nucleon. 

The effective field operator V(aS) has a pole in the ex- 
citation energy a) s of the state \s) under consideration, 



(v A(a) s )g 0s ) g 0s 
V(u) = X - '- + V r (oj). 
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Fig. 1. The proton and neutron normal transition amplitudes g {0) 
in " 8 Sn nucleus. 
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Fig. 2. The charge transition densities p^ 0) in " 8 Sn nucleus. Red 
line corresponds to surface pairing, blue one, to volume pairing. 
Experimental data are taken from 1131 . 



The quantity go s has the meaning of the corresponding ex- 
citation amplitude. It obeys the homogeneous counterpart 
of Eq. (0 and is normalized as follows J4), 



g Qs ^\ = -i. 



(ii) 



with obvious notation. 

It should be noted that the above TFFS equations are 
essentially equivalent to QRPA those. Important difference 
occur if one goes beyond the QRPA IfTTI . For the pairing, 
we use the diagonal approximation, A^=A^6u', and "de- 
veloped pairing" formalism which neglects effects of par- 
ticle non-conservation. According to 0121 . it leads to errors 
^0.1 MeV in values of Ax- 



3 Surface versus volume pairing for 
semi-magic nuclei 

The QRPA-like equations of the self-consistent TFFS were 
solved in the coordinate representation, with exact consid- 
eration of continuum. We limit ourselves with semi-magic 



nuclei where the phonon coupling corrections to QRPA are 
less than in nuclei with both non-magic subsystems. As a 
benchmark nucleus, we take I18 Sn which is in the middle 
of the tin chain. Let us begin from the normal components. 
As it is seen from Fig.[T] the normal transition amplitudes 
g(°\r) found for these two kinds of pairing are rather close 
to each other. The same is true for the normal components 
of transition densities that is illustrated for the charge tran- 
sition density p ^ (r) which is compared in Fig. [2] with 
the experimental one found in Ref. Ifl3l with the model- 
independent analysis of the high-precision data on the elas- 
tic electron scattering. It is worth mentioning that these two 
figures confirm the interpretation by Khodel [14. 8 1 of the 
2|-excitations as quantum capillary waves belonging to the 
Goldstone branch of surface vibrations. 

Another situation appears for the anomalous transition 
amplitudes gr (1 ' 2) , see Fig. [3] We see that now these compo- 
nents in the case of the surface pairing are approximately 
3 times stronger than for volume pairing. Unfortunately, 
there is no direct way to observe the anomalous transi- 
tion amplitudes </ 1,2) , but there is an indirect one. When 
comparing Fig.Q]and[3] one should take into account that 
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Fig. 3. The neutron anomalous transition amplitudes g {l2) in 
118 Sn nucleus. Red solid and dotted lines correspond to surface 
pairing, blue ones, to volume pairing. 



the amplitudes gr (0) and gr (1 ' 2) have the same normalization, 
therefore their absolute values may be compared directly. 
Thus, we see that in the nucleus under consideration the 
anomalous components are comparatively small, approx- 
imately 3 times smaller than the normal one even for the 
surface pairing case. Therefore we can consider the anoma- 
lous components as a small perturbation and present the 
eigen-energy of the L-phonon as 



(°) , x (!) , s; P) 

(jl>i-u l +otir L +ocl> l , 



(12) 



where the corrections due to the anomalous components 
are 

^ *-(«/«>)>,. (13) 

Here gr (1 ' 2) denote average values of the corresponding ma- 
trix elements. Therefore, for the surface pairing this (neg- 
ative) correction should be larger. It turns out that way in 
actual calculations. For B(E2) values, there is an interfer- 
ence between the normal amplitude and the anomalous 
ones, g (1 ' 2) , therefore dependence of the result on the type 
of pairing is not so definite. 

Calculation results for a>2 values in tin isotopes are dis- 
played in Fig. [4] We see that for the surface pairing a>2 val- 
ues are approximately by 0.3 Me V less than for the volume 
one, in a qualitative accordance with Eq. (fT3l . The surface 
version turns out to be closer to the experiment. The rms 
deviation of the theory from the experiment is Au> rms =0.16 
MeV for the surface pairing and /Jw lms =0.37 MeV for the 
volume pairing. When calculating the rms values of Acj^, 
we excluded the double magic nuclei and their neighbors, 
1()2 Sn and 130_134 Sn. Magic nuclei are excluded on the rea- 
son of absence of pairing and hence the effect under con- 
sideration. As to their neighbors, the reason is deficien- 
cies of the "developed pairing" approximation we use for 
such nuclei. First, the particle number conservation only 
in average inherent to this approximation works for them 
rather poorly, as it follows from the analysis in Ref. [ 12 1. 
Second, it is true for the equality of the operators A + and 
A ~, see Ref. ||4), which is supposed in this approximate 
scheme. Indeed, e.g. for an isotope with the neutron num- 
ber N = N m - dg - 2, the operator A + connects the nucleus un- 
der consideration with the one in which the neutron pairing 
is absent. 
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Fig. 4. Excitation energies w(2|) for tin isotopes. Experimental 
data are taken from 1151 . 

The B(E2) values in the tin chain are shown in Fig. [5] 
In this case, the situation is not so definite because of ex- 
perimental errors. It is worth mentioning that in this case 
the volume version looked a bit better when the old data 
were used J2]], but new data of Ifl6ll agree better with the 
surface version. A comparison is made with the Skyrme- 
Hartree-Fock-Bogolyubov (SHFB) theory predictions of 
Ref. 03 with the SkM* and SLy4 Skyrme force. Among 
these two SHFB results, the SkM* ones are significantly 
better. Our results for the excitation energy, especially for 
the surface pairing, are much better than both the SHFB 
ones. The same is true for the B(E2) values. A general de- 
ficiency of the two SHFB theory predictions is their non- 
regular A -dependence whereas the data depend on A rather 
smoothly. 

For the lead isotopes, the situation is qualitatively sim- 
ilar, see Figs. [6] [7] Again, a>2 values for the volume pair- 
ing are at 0.3 MeV higher of those for the surface case. 
Now we have zlw rms =0.33 MeV for surface pairing and 
Aa> rmti =0.47 MeV for volume pairing. In accordance with 
the above arguments, the calculation of the A(L> rmti values 
is carried out with exclusion of 206 - 210 pb nuclei. We see 
that here accuracy of our calculations is worse than for the 
tin chain but again the surface pairing model predictions 
look better. Here the SHFB results for both the versions 
of Skyrme force look much better than for the tin chain 
and their accuracy is approximately the same as ours for 
the surface pairing case. The B(E2) values are known only 
for four lead isotopes, 204 - 210 Pb, the SHFB accuracy is a 
bit higher. Note that the 204 Pb nucleus is the only among 
them which is described within our QRPA-like scheme ad- 
equately. 

Let us go to the nickel isotopes. Fig. [8] demonstrates 
that the two versions of the pairing force we compare are 
practically equivalent to each other from the point of view 
of reproducing the data on the neutron separation energy 
S„ for the nickel chain. Analogous demonstration for the 
tin and lead isotopes can be found in [2|. Calculation re- 
sults for excitations energies are displayed in Fig.|9]and for 
B(E2) values, in Fig. [10] The chain contains two double- 
magic nuclei, 56 Ni, N = 28, and 78 Ni, N = 50. In addition, 
the nucleus 68 Ni with the neutron number N - 4-0 also 
can be interpreted as a quasi-magic one. Indeed, the exci- 
tation energy of the 2\ state in this nucleus is significantly 
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Fig. 5. 6(£2,up) values for tin isotopes. The "old" experimental data are taken for 114 ~' 24 Sn from |T5), for 126 ~ 134 Sn from (T8|, and 
for 106 ~ U2 Sn from [191201211 . The new values marked "GSI" and "ANU" are taken from 1 16|. The first values are measured at the at 
the UNILAC accelerator of the Gesellschaft fiir Schwerionenforschung (GSI). The second ones are found also in 1 16] by analyzing data 
obtained in the magnetic moment measurement 1221 at the Australian National University (ANU). 
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Fig. 6. Excitation energies w(2j) for lead isotopes. 
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Fig. 7. B(E2, up) values for lead isotopes. 



higher than in neighboring even nickel isotopes that is a 
typical feature of magic nuclei. Although the effect is not 
so strong as for 56 Ni, it is also rather pronounced. For a 
complete magicity, the distance between occupied single- 
particle levels and the next free one should be greater than 
the double average gap value 2A. In this case, the gap equa- 
tion has only the trivial solution. In the case under discus- 



sion, the distance between the occupied 2p - l/-shell and 
the next level lgg/2 is of the order of 2/1 and the gap equa- 
tion has a non-trivial solution. However, some features of 
magic nuclei occur in this nucleus but in a smoothed form. 
Thus, just as for the tin and lead chain, for analyzing the 
"surface versus volume" effect under discussion, it is worth 
to exclude the neighboring to 56 Ni and 78 Ni nuclei. In ad- 
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Fig. 8. Neutron separation energies S„ for Ni isotopes. 




Fig. 9. Excitation energies &>(2t) for Ni isotopes. 
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Fig. 10. B(E2, up) values for Ni isotopes 

dition, the isotopes 66 ~ 70 Ni should be also considered with 
some caution. 

We see that qualitatively the pattern is similar to that 
in the tin and lead chains, however the effect for a>2 value 
is now not so regular. Indeed, for the isotopes 72 ' 74 Ni it is 
practically zero and for the isotope 7() Ni it is of the oppo- 
site sign, although rather small. Although the latter could 
be interpreted as the effect of closeness to the quasi-magic 
nucleus 68 Ni, let us try to understand the reason of such 
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Fig. 11. The proton and neutron normal transition amplitude g (0) 
in the 70 Ni isotope. 
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Fig. 12. The neutron anomalous transition amplitudes and 
g {2) in the 70 Ni isotope. Red solid and dotted lines correspond to 
surface pairing, blue ones, to volume pairing. 



contradiction to the formula (12) within the scheme we 
use. For this aim, we displayed for this nucleus the nor- 
mal and anomalous transition amplitudes in Fig. QT| and 
Fig- 02] correspondingly. We see that, just as for the I18 Sn 
nucleus, the normal amplitudes for the two models of pair- 
ing under consideration are very close to each other. In 
addition, again the anomalous amplitudes for the surface 
pairing case are much stronger than those for the volume 
pairing. However, in this nucleus the absolute values of the 
anomalous amplitudes are rather close to those of the nor- 
mal amplitudes. Therefore, the consideration of the anoma- 
lous amplitudes as a small perturbation, which is the base 
for receiving the estimation (12), is not valid now. 

Agreement with the data on 102 values for the surface 
pairing case is almost perfect for nuclei lighter than 68 Ni, 
however it becomes worse for heavier isotopes. The most 
strong disagreement for 7() Ni and 76 Ni nuclei again can 
be interpreted as the effect of closeness to magic nuclei 
(a quasi-magic in the first case). One more reason exists 
for poorer agreement for heavy nickel isotopes. For them, 
neutrons are filling up an almost isolated level 1^9/2. In- 
deed, as it was already discussed, the distance to the filled 
levels of negative parity is about 2A whereas that to neigh- 
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boring levels of positive parity is significantly greater than 
2A. In such a situation, predictions for the positive par- 
ity states depend significantly on details of the calculation 
scheme. The experimental values of transition probabili- 
ties are known only for isotopes till 68 Ni, here agreement 
is quite reasonable. 



4 Conclusion 

Excitation energies, transition probabilities and other char- 
acteristics of the first 2 + excitations in three chains of semi- 
magic isotopes, the lead, tin and nickel ones, are calculated 
within the self-consistent TFFS on the base of the EDF 
by Fayans et al. |6). The DF3-a functional IflOl is used 
which differs from the original one DF3 with the spin-orbit 
and effective tensor terms only. A reasonable agreement 
with available experimental data is obtained. Predictions 
for other isotopes are made. 

The effect of the density dependence of the effective 
pairing interaction is analyzed by comparing predictions 
of two models of the effective pairing interaction, the sur- 
face pairing and the volume one. The effect is found to be 
noticeable, especially for the 2 + -energies which are, with 
the only exception, systematically higher at 200-300 keV 
for the volume paring as compared with the surface pair- 
ing case. Thereby predictions of the surface pairing model 
agree with the data better than those of the volume model. 
Thus, the analysis of the excitation energies of the first 2 + - 
states in semi-magic nuclei gives a new evidence in favor 
of the surface pairing. As to B(E2) values in tin isotopes, 
the volume pairing model looked better in comparison with 
old experimental data, but the new data of lfl6l agree better 
with the surface model. 
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